Time-Resolved Upconversion of Entangled Photon Pairs 
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In the process of spontaneous parametric downconversion, photons from a pump field are converted 
to signal and idler photon pairs in a nonlinear crystal. The reversed process, or upconversion of 
these pairs back to single photons in a second crystal, is also possible. Here, we present experimental 
measurements of the upconversion rate with a controlled time delay introduced between the signal 
and idler photons. As a function of delay, this rate presents a full width at half maximum of 27.9 fs 
under our experimental conditions, and we further demonstrate that group delay dispersion of the 
photon pairs broadens this width. These observations are in close agreement with our calculations, 
thus demonstrating an ultrafast, nonclassical correlation between the signal and idler waves. 

PACS numbers: 42.50.Dv, 03.65.Ud, 42.65.Lm, 42.65.Ro 



There are few ways of characterizing photon wavepack- 
ets on an ultrafast time scale. While temporal resolution 
of a few fs can be necessary, optical detectors have a re- 
sponse that is slower by, at least, several orders of magni- 
tude. One way of overcoming this limitation is with two- 
photon interference, as pointed out by Hong, Ou, and 
MandelQ. With fs-level resolution, they determined the 
distribution of arrival time intervals between the signal 
and idler photons produced by spontaneous parametric 
downconversion (SPDC). Refinements of this approach 
have since measured relative signal/idler group delays to 
a remarkable precision of 0.1 fs [2]. 

Other interactions produced with nonclassical light 
have been considered in previous works [1, 0, 0, @, 13|- 
These interactions inspire another possible ultrafast char- 
acterization method: to introduce a controlled delay be- 
tween the signal and idler photons of SPDC, overlap 
them spatially in a nonlinear crystal, and detect the sum- 
frequency photon produced should they upconvert. This 
approach would yield high-resolution temporal informa- 
tion about the original quanta; indeed, it is analogous 
to the widely-utilized autocorrelation method of char- 
acterizing classical laser pulses Q, where fs-level resolu- 
tion is common. Moreover, it will be shown here that 
such an approach is dispersion-sensitive, in contrast to 
the Hong-Ou-Mandel interferometer Some years ago, 
this type of time-resolved upconversion was studied with 
high-gain (w 10^°) parametric light However, it was 
soon pointed out that the approach would be of no use 
at the photon level[l[, essentially due to the weak non- 
linearities of crystals then available. 

Present-day crystals have stronger nonlinearities and 
it is therefore time to readdress this issue. Using an un- 
usually bright SPDC source, Dayan et have pro- 
vided the first observations of signal/idler paired upcon- 
version in a nonlinear crystal, where photon pair mem- 
bers were kept perfectly synchronized. Yet, two recent 
theoretical works have considered time-resolved, pairwise 



upconversion [d. Ml , and related experimental results have 
been reported[lO| that we will discuss later. In part, our 
purpose is to present novel, time-resolved experimental 
results, but we begin with theoretical considerations. 
The two-photon state of SPDC may be written as 



|*Dc) = \0)+V 



d^hfociks, h)al{ks)al{h)\0), 

where ry is related to the conversion efficiency, fociks, ki) 
is the joint amplitude, a^^{k^) is the creation operator for 
the signal (/i = s) or idler (n =i) modes with wavevector 
kfj,, and |0) is the vacuum state. We want to describe the 
upconversion of these photon pairs into single photons at 
the sum frequency. In the spontaneous conversion limit, 
the upconverted state \'^uc) can be expressed in the form 



|*C/C) = (l + C/l)|*Dc) = I 1 + ^ 



dt't'Hit') I \^DC), 
(2) 

where (1 + Ui) is the unitary time-evolution operator. 
The upconversion Hamiltonian H(t) can be expressed in 
terms of the nonlinearity d{f) and electric field operators 

El'^\r,t) as 

Hit) = I d'r d{r) (f, t) Ei+^ (f, t) (f, t) + H.C. , 

(3) 

where fi = u now also includes the upconverted mode, 
-I-/- denotes positive-/negative- frequency parts, and the 
second term is the hermitian conjugate. By expressing 
each of the electric fields as a Fourier integral in k^, and 
performing the integral of Eq. ([2]), we obtain 

Ui ^ d^kuj d%J d%f^ciks,kku) 

X S{u!u - ^^s - i^i)dl{ku) dsiks) diih) . (4) 
Defined in terms of fc-|_^ = ksu + k^i, (with v ~ x,y), 
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fuc{ks, ki, ku) follows from the integral of Eq. [3]as 

fuc{ks,ki, ku) = 5{k+.j,-kux)5{k+y~kuy)<^uc{ks,ki, ku) 

where <^uc{ks,ki,ku) is given by sinc(/3) cxp(— with 
(3 = /\kzL/2 and AA:^ kp^ - k^z - - kg, L 
is the upconversion crystal length, and kg ~ 2n/A 
is the wavenumber associated with the crystal poling 
period A. We now assume that the SPDC pump is 
monochromatic with frequency ujp. In this case we 
can write fDc{ks,ki) as 5{ijJs+uJi—ijjp) g{ks,ki) and the 
upconverted state of Eq. ([2|) becomes \^uc) = |0} + 
K / d^ku F{ku)al{ku ,ijJp)\Q) . This integral is over only 
transverse components of km the creation operator la- 
bel is now expressed in terms of k^ and LOp, k includes 
various constants, and i^(fc^) is given by 



F{k^) 



(f fcs $^(7 i^s ,k^;ujp - ujs, kl ; ujp , fc,„ ) 



X 9{^^s-,k^;ujp- ujs-,k^) , (6) 



where both ^uc{') a-nd g{-) are now written in terms of 
the new variables. 

For a plane- wave SPDC pump, g{LOs, k;^; uip — uis, ku) is 
given as d{k+x) S{k+y) ^ociks, h) where 3>_Dc(fcs,fcj) = 
sinc(/3) exp(i/3), and we have assumed that both crystals 
are characterized by the same /3. We further assume that, 
before the upconversion crystal, the signal and idler pho- 
tons propagate through dispersive optical elements giving 
each a spectral phase exp[i(/)^(a;^)], that they encounter a 
limiting aperture described by pupil function P{k^), and 
that the signal mode is delayed by t (these operations 
amount to multiplying fDc{ks,ki) of Eq. ([1]) by appro- 
priate factors, which simply carry through the analysis 
above). The upconverted state may then be written as 
|*r/c) = |0) + K'|«'g^) with 



I* 



uci 



duJsS{LUs)e 



at(0,wp)|0), (7) 



where 



d2fc^J-sinc2(/3)e*^'^'-^e'['^=("=)+'^-("-'lp(fc,)P(fc,), 



(8) 

and overlined quantities are evaluated with = o^p — lUs 
and k^ = —k^. 

Our main theoretical results are Eqs. ([TJE])- Here, 
S{ijJs) is related to the amplitude spectrum of signal- 
mode photons that participate in the upconversion pro- 
cess. Still, it is also clear that S{ijJs) depends on the 
dispersion imposed on both the signal and idler modes 
through the second exponential factor of Eq. ([8]). At the 
least, this indicates that control of dispersion is necessary 
in the experimental work to be described. We also note 
that previous works 0| have derived what amount to 
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FIG. 1: Experimental apparatus. Parametric downconver- 
sion from lithium niobate crystal LNl is coUimated and sent 
through a series of 4 Brewster prisms, to be upconverted in 
a similar crystal LN2. Two mirrors (one fixed, the other on 
a piezoelectric stage, translated as shown) introduce a time 
diff'erence r between the upper and lower optical paths. 



axial limits of Eqs. (I7][8]), and for us it is also straight- 
forward to obtain a simplification of Eq. ([5]) in the axial 
case. However, the factors related to phasematching are 
strongly angle-dependent in cases later considered here 
and, to obtain accurate results, retaining the form of 
Eq. ([8]) is necessary. 

The experiment is shown in Fig. 1. The pump laser 
was a single-mode, continuous-wave ring laser of wave- 
length 532 nm and power 1 W. Its beam was focused to 
a 45 /im waist in an uncoated, magnesium oxide doped, 
5 mm long lithium niobate crystal, which produced the 
down-converted light. The crystal was periodically-poled 
and temperature-controlled to permit phasematching to 
nearly axial, co-polarized, frequency-degenerate photon 
pairs. The emission within a cone of half-angle 2° was 
accepted by the aperture of the optical system that fol- 
lowed. 

The light was coUimated by a 75 mm focal length lens 
and was sent through a SF14 Brewster prism, which de- 
viated the pump beam from the system. The down- 
converted light passed through four such prisms, sym- 
metrically arranged, to compensate for the chirp of the 
photon wavepackets arising from all dispersive media in 
their pathQ. The tip-to-tip spacing of the first or second 
prism pair was 352 mm. 

In the coUimated downconverted beam, transverse mo- 
mentum conservation implies that a given photon will 
have its pair member on the opposite side of beam cen- 
ter. Thus, in the light exiting the prisms of Fig. 1, re- 
flecting the upper part of the beam (containing, say, sig- 
nal photons) from a displaceable mirror and reflecting 
the lower part (idler photons) from a second, fixed mir- 
ror provided the signal/idler time difference r. At beam 
center, a 1.5 mm intermirror gap guaranteed that any sig- 
nal or idler photon could interact with only one mirror, 
thus dividing reflected pair members into two distinct, 
time-shifted modes. Angular mirror co-alignment was 
achieved with a subsidiary Michelson interferometer (not 
shown), where the mirrors acted as a single end-mirror. 
There, a tunable source was used to find the translation 
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FIG. 2: Experimental results, (a) Upconversion rate R{t) 
versus signal/idler delay r for optimum dispersion, (b) R{t) 
with upconversion crystal 20° C above phasematching tem- 
perature. Other results are for additional signal/idler paths 
in 6mm of fused silica (c, solid curve), 12mm of fused silica 
(c, dot-dashed curve), and 37mm of SFIO glass (d), show- 
ing increasing broadening of R{t). Data in (c)-(d) have been 
scaled to correct for Fresnel window losses. 

setting that made the mirrors coplanar (t = 0), for which 
fringes remained ahgncd across the gap during tuning. 

The light was then focused in the second crystal, which 
had antireflection coatings but was otherwise identical to 
the first. The upconverted flux passed through a BG39 
green filter, was focused into a multimode fiber, and 
was sent to a SPCM-AQR-13-FC single-photon count- 
ing module. Using prism insertion to adjust dispersion, 
the upconversion rate for r = was maximized (hereafter 
referred to as 'optimum dispersion'). This maximum 
rate was found to occur for 1.5° C less than the axial 
degenerate-phasematching temperatures of the crystals 
(w50° C), which slightly opens the phase-matched emis- 
sion and acceptance cones. In our results, data are the 
count rate obtained over 6 s and, at the close of each data 
set, the background rate (detector dark count plus stray 
light, «195s~^) was measured over 30s and subtracted. 

In Fig. 2(a) we show results for the upconversion rate 
R{t) with optimum dispersion, where a rate of nearly 
800 s~^ appears at r = 0. This peak decays with width 
(full width at half maximum throughout) 27.9 fs, and 
minima and secondary maxima appear near, respectively, 
r = ±30fs and ±44fs. For |t| >70fs, R{t) decays to a 
noise floor of essentially zero mean. While statistical er- 
ror bars are too narrow to be shown clearly, they are 
consistent with the residual fluctuations seen (error bars 
of ±13s~^ near r = 0, and ±6s~^ in the noise floor, 
for example). This low noise floor indicates that random 
upconversion (i.e., that between unpaired photons) is not 
significant under the experimental conditions. As a check 
on the validity of these data. Fig. 2(b) shows data taken 
under identical conditions, but with phasematching sup- 
pressed in the upconversion crystal by raising its temper- 
ature by 20° C. Here R{t) has fallen to a noise level of 
essentially zero mean, further supporting that the pho- 
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FIG. 3: Theoretical results, (a) Upconversion rate R{t) (ar- 
bitrary units) versus signal/idler delay r for optimum dis- 
persion. Also shown is R{t) for additional common sig- 
nal/idler dispersion due to 6 mm of fused silica (b, solid 
curve, <j!." = 99fs2), 12 mm of fused silica (b, dot-dashed curve, 
(?!>" = 198 fs^), 5mm of SFIO glass (c, 0" = 513 fs^), and 37mm 
of SFIO glass (d, <?!>" = 3790 fs^). 



tons of Fig. 2(a) arose from the desired, phase-matched 
process. 

The consequences of the temporal stretching of the 
photon wavepackets are shown in Figs. 2(c)-(d). Return- 
ing to the optimum case of Fig. 2(a), windows of fused 
silica were then inserted just before the fourth prism of 
Fig. 1, spanning the signal and idler paths, to intro- 
duce an identical group delay dispersion 0" in each mode 
(taken as the second derivative of accumulated window 
phase, calculated at degenerate frequency). In Fig. 2(c) 
(6 mm thickness of fused silica, (jj" = 99fs^), the peak of 
R{t) has fallen to w 640 s~^ and has widened slightly to 
29.4 fs. Using a similar window of twice this thickness 
(0" = 198 fs^), also in Fig. 2(c), the upconversion peak 
continues to fall (Ri 420s~^) and widen (36.0 fs width). 
Finally, we have the extreme case of Fig. 2(d), where a 
window of 37 mm of SFIO glass {(j)" = 3790 fs^) was em- 
ployed. A low, featureless result is obtained at a level of 
20 ± 6s~^, indicating that R{t) has been so broadened 
that pairwise upconversion is scarcely observable. 

We now compare with theoretical results. Using the 
experimental parameters, R{t) was computed from the 
squared modulus of the term in the square brackets of 
Eq. [7] with 4>s{uj) = '/'i(w), the integral being done nu- 
merically, and the optical system dispersion being cal- 
culated exactly. As in the experiment, prism insertion 
was varied to find optimum dispersion (i.e., maximum 
i?(0)); this occurred for a residual system (from first- to 
second-crystal centers) group delay dispersion of 28 fs^, 
which compensates for quartic dispersion. In Fig. 3(a), 
this optimal R{t) presents a main peak of width 25.0 fs, 
with secondary maxima at T = ±42fs. Fig. 3(b) shows 
R{t) for the additional signal/idler dispersion of 6 and 
12 mm of fused silica, which lowers the peaks and broad- 
ens the widths to, respectively, 26.0 fs and 30.9 fs. This 
trend continues with increasing dispersion in Figs. 3(c) 
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and 3(d) to the extent that peaks no longer appear at 
r = 0. The comparison with our experimental results is 
quite good, as may be seen by comparing Figs. 3(a), 3(b), 
and 3(d) with plots at similar spatial positions in Fig. 2. 
The most significant difference is that the calculated R{t) 
is slightly narrower; much of this difference arises because 
the theoretical single-photon SPDC bandwidth (130 nm) 
is broader than that experimentally measured (115 nm), 
which implies shorter wavepackets in the calculations. 

Experimental results that parallel ours, but that have 
significant differences, have been reported elsewhere (lo|. 
There, a V-shaped spectral phase mask centered at the 
degenerate frequency was applied to both photons, pro- 
ducing a relative time shift between the lower- and 
higher- frequency pair members. It is readily shown that 
this amounts to replacing the factor exp(— iwgr) in Eq.[7| 
with exp(— i \llIs — Wp/2| r) which, for a given S{ujs), pro- 
duces a different upconversion trace R(t). For example, 
with a Gaussian S{ujs) our Eq. W\ predicts a Gaussian 
R{t), while the approach of Ref. Il0| predicts a curve 1.7 
times wider, with tails similar to a Lorentzian. Moreover, 
it is readily shown that the approach of Ref. [l^ predicts 
skewncss in curves analogous to Fig. 3(b), since disper- 
sion produces time shifts between the lower- and higher- 
frequency pair members (although we do not reproduce 
these calculations here). In contrast, our experiment in- 
stead relies on spatial separation to define the signal and 
idler modes, and the spectral content of each mode is left 
intact. Thus, while our method produces a pure time 
shift between the two modes, the approach of Ref. [13 
implies a temporal re-shaping of the single photons. 

It is also notable that upconversion of single photons 
with a classical laser beam has been demonstrated [ll| . 
and an extension of this approach has been reported 
in which each member of a given signal/idler pair up- 
converts with the same ultrafast laser pulse [l3l. What 
is obtained, after averaging, is the two-dimensional his- 
togram of signal and idler arrival times with respect 
to the laser pulses, with resolution determined by the 
pulse width (150fs)[l3|. In comparison, our method 
is intrinsically only one-dimensional in time, although 
it has at least 2 orders of magnitude more resolution, 
depending simply on mirror positioning accuracy. Fi- 
nally, it is significant that our approach is highly sen- 
sitive to group delay dispersion of the photon pairs 
while, in the similar case of a monochromatic SPDC 
pump, the Hong-Ou-Mandcl interferometer is not[9|. In 
fact, in the phase term of Eq. ([5]), the phase factor 
exp {iAkzL/2) exp {i[(ps{ij->s) + is that associated 

with the dispersed two-photon state itself, and thus our 
approach provides a fairly direct, high-sensitivity means 
of probing the consequences of its dispersion. Due to 
slow detector response, observing group delay dispersion 
effects through direct detection of signal and idler arrival 
times has required « 10^ times more dispersion than we 
easily see, for example, in the cases of Figs. 2(c) and 



3(b) [l3|. Based on the sensitivity of our calculations, we 
estimate that wlOfs^ of group delay dispersion could be 
observed with our method (in either the signal or idler 
mode, or distributed in both, as is clear from the phase 
terms of Eq. ([5])). However, an experimental demonstra- 
tion of this sensitivity is beyond the scope of our work. 

In conclusion, we have reported results in which entan- 
gled photon pairs are upconverted with a controlled time 
delay introduced between the signal and idler photons. 
We have thus demonstrated an ultrafast, nonclassical 
correlation between the signal and idler waves, even for 
our case of a continuous- wave, single-mode SPDC pump. 
Our theoretical analysis has developed an expression for 
the upconverted amplitude that is strongly dependent on 
the dispersion of the signal and idler modes. In experi- 
ments, the signal photons have been spatially separated 
from the idlers, delayed with the fs-level resolution inher- 
ent in a piezoelectric stage, and then recombined with 
the idlers for upconversion. Close agreement has been 
demonstrated between the observed and predicted up- 
conversion rates as a function of signal/idler time delay. 
We believe that this technique is a valuable tool for the 
characterization of entangled photon pairs at the fs level. 

This research was supported in part by CONACYT 
(Mexico) under Grant 49570-F. 
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